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I . Objectives 

The objectives of research performed under this grant 
are (1) to develop techniques to detect and remove long-term 
instrument drifts from the operational SBUV/2 data base 
using periodic measurements made by the Shuttle Solar 
Backscatter Ultraviolet Spectral Radiometer (SSBUV) and (2) 
to evaluate the potential effect of scattering by 
atmospheric aerosols on our ability to detect trends in 
ozone using SBUV-class instruments. Work during the first 
six months of this grant, beginning in mid-February 19B7, 
has focused on the first of these objectives. 

I I . Methodology 

A method to combine backscatter albedo measurements 
from SSBUV and the SBUV/2 instruments for the purpose of 
removing drifts in the latter data set has been developed 
and tested. At present the techniques have been applied to 
a simulated data base which represents the two year lifetime 
of a single SBUV/2 instrument. The following paragraphs 
summarize the mathematical model. 

The backscatter albedos, A(i) = I(i)/F(i), are the 
basis for SSBUV-SBUV/2 comparisons, where I < i ) is the 
backscattered radiance which emerges from the atmosphere in 
the vertical direction, F(i) is the incident solar 
irradiance, and the index " i " labels wavelength. During a 
single flight of SSBUV, denoted by the index " j " , there will 
be k=l,2,...,K coincident measurements with the operational 
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SBUV/2. Let A„(i,j,k) denote the albedos measured by SSBUV 
and A = (i,j,k) refer to nearly simultaneous SBUV/S data. 
During a given two year SBUV/2 flight j will take on values 
j=l,2,3, and perhaps b . A typical value for K, the number 
of SSBUV-SBUV/2 coincidences on a single Shuttle flight, is 
K=30. 

We define a "correction factor" c(i,j) which normalizes 
the SBUV/2 albedos to those measured by SSBUV. This is: 

K 

c(i,j) = E CA.(i, j,k)/A e (i, j,k)l/K (1) 

k — i 


The statistical uncertainty in the correction factor is 
defined by the variance: 

K 

<r e = E CA.(i, j,k)/A=(i, j,k)-c(i, j)D s /(K-l) (2) 

w — i 


The SSBUV-SBUV/2 intercomparisons made on a given day define 
one correction factor for each wavelength (index " i " ) . This 
could be altered in the future if, for example, one wished 
to derive a single correction factor for an entire SSBUV 
mission lasting several days. Given several flights of 
SSBUV during the lifetime of a single SBUV/2 instrument, one 
can derive a time dependent correction factor by fitting a 
linear or quadratic curve to the c(i,j), j=l,2,3, and maybe 

*t. 


A 


or 


c(i,t) = c 0 <i> + Ci<i)t 


c(i»t) = c 0 (i) +■ Ci<i)t +c e (i)t E 


(3a) 


(3b) 


where "t" is time in days from the start of data collection 
by the SBUV/2. A "corrected" set of SBUV/2 albedos, A c , as 
a function of time is then: 


A c (i,t)=c(i,t)A e (i»t) (^) 

These corrected SBUV/2 albedos are the data set for use in 
seeking long-term trends in atmospheric ozone. 

A central issue of this study involves definition of 
how uncertainties in the measured albedos and a lack of 
exact simultaneity between SSBUV and SBUV/2 influence the 
derived correction factors. To investigate this we have 
developed simulated data sets for A m and A = which contain 
known systematic and random errors. 

The Nimbus 7 SBUV instrument has provided a set of 
measured albedos which, for test purposes, we take to be the 
true atmospheric albedo with zero error. We consider data 
measured over the latitude range 0 to 15°N here. We then 
generate a two year series of simulated SBUV/2 albedos, A E > 
via : 


A s ( i , t ) 


Cl+f ( i ,t )+€( t)]A T ( i ,t) 


(5) 
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Here A T (i,t) is based on the monthly mean albedo from Nimbus 
7 plus a random noise component. Daily values of A T <i»t) 
are chosen at random from a normal distribution whose mean 
is the average Nimbus 7 albedo for the month and whose 
standard deviation (<r) is 0.333*/. of this average (3 a = 1*/.). 
The function f(i,t) simulates long-term drifts in the SBUV/2 
instrument. We have examined both linear and exponential 
drifts: 


f(i,t)=d(i)t (6a) 

and : 

f ( i » t ) = d ( i ) { 1 - expC-t/T ( i ) ] > (6b) 

where we specify d(i) and T(i). Equation 6b simulates an 
instrument which experiences an initial drift but which 
eventually stabilizes over a time t equal to several times 
T(i). The parameter €(t) provides for noise in the 
measurements. We select each value of € at random from a 
normal distribution whose mean is zero and whose standard 
deviation is a = 0. 004/3. 

We also simulate a set of SSBUV albedo measurements by: 
A.(i,j,k> = C 1 +b ( j ) 3 Cl+e~(k)3 A T ~(i,j,k) (7) 

The index " j " labels a flight of SSBUV, and "k" labels a set 


of albedo measurements obtained nearly simultaneous with the 
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SBUV/2. Taken together j and k specify a unique time t in 
equation 5 for which a comparison between SSBUV and SBUV/2 
is possible. The albedo A T ~ in equation 7 is selected from 
a normal distribution whose mean is the Nimbus 7 SBUV albedo 
and whose standard deviation is 0.333V.. Note that 
A T *<i,j,k) in equation 7 is not numerically identical to 
A-r<i»j»k) in equation 5. The common point is that both A T 
and A T * are selected from the same normal distribution based 
on Nimbus 7 data. The random selection of A T and A T ~ from 
the same distribution simulates a lack of exact simultaneity 
between SSBUV and SBUV/2 and the fact that the fields of 
view of the two instruments cover a different area when 
projected onto the earth’s surface. 

The term €~<k) accounts for SSBUV repeatability and is 
determined in the same manner as C(k) for SBUV/2. The 
parameter b(j) allows for calibration biases in SSBUV which 
can vary from one flight to the next. The sample 
calculations presented in this report select b(j) at random 
from a normal distribution with a mean of zero and a 
standard deviation of 0.333V. <3 ct = IV.). Future calculations 
will specify a known b(j) for each Shuttle flight to examine 
the sensitivity of the correction factors to SSBUV 
calibration biases. 

III. Sample Results 

Figures 1 through A present samples of our preliminary 
results. The top panel of each figure gives the simulated 
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SBUV/2 albedo data set over a two year period. The origins 
of the time scales (expressed in days) refer to the start of 
SBUV/2 data collection. The jagged behavior in the 
simulated SBUV/2 data set arises from the use of monthly 
mean albedos as opposed to daily values. This has no impact 
on the derived correction factors. The stars in the top 
panels indicate SSBUV measurements. In Figures 1 and 2 we 
assume three SSBUV missions separated by eight months, while 
Figures 3 and ^ contain four SSBUV flights separated by 
seven months. Each star is the mean of six albedos obtained 
at the SSBUV-SBUV/2 coincidences on a single day, and we 
assume a five day Shuttle mission. The stars in the center 
panels of all figures represent correction factors derived 
for each day of a SSBUV mission. 

The solid line is a fit of equation 3a (Figures 1 and 
2) or 3b (Figures 3 and **) to the individual c(i,j) values. 
The bottom panel of each figure compares the true SBUV/2 
drift, Aj=— A t (darker line), with the computed drift, A C -A T 
(lighter, dotted line). If the correction procedure were 
perfect the two lines would coincide. Figures 1 and 2 refer 
to wavelengths 273.6 and 2B7.7 nm respectively and assume a 
linear SBUV/2 degradation over time at a rate of 4*/. per year 
(d(i) = O.O^t/365 in equation 6a). Figures 3 and 4 refer to 
an exponential instrument drift followed by stabilization, 
as described by equation 6b with d(i) = 0.0A and T(i) = 182 
days. These results should be regarded as preliminary, 
although they indicate that SSBUV is capable of detecting 


and correcting SBUV/2 drifts of the magnitude assumed here 
to a high degree of accuracy* although the procedure is not 
totally free of error. Many more test cases must be 
analyzed here before definitive conclusions are reached. 
This activity will continue in the coming six month period. 
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Figure 3 


Wavelength: 273.61 
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